Controllable Gaussian-qubit interface for extremal quantum state engineering 
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We study state engineering through bilinear interactions between two remote qubits and two-mode Gaussian 
light fields. The attainable two-qubit states span the entire physically allowed region in the entanglement- 
versus-global-purity plane. Two-mode Gaussian states with maximal entanglement at fixed global and marginal 
entropies produce maximally entangled two-qubit states in the corresponding entropic diagram. We show that a 
small set of parameters characterizing extremally entangled two-mode Gaussian states is sufficient to control the 
engineering of extremally entangled two-qubit states, which can be realized in realistic matter-light scenarios. 
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The structure of quantum correlations within a given sys- 
tem depends strongly on the dimension of the state spaces of 
its constituents. The relation between correlations and system 
dimensionality becomes particularly relevant when consider- 
ing states of bipartite systems whose parties are defined on 
Hilbert spaces of different dimension. An extreme instance is 
that of a bipartite compound consisting of a continuous vari- 
able (CV) system with infinite-dimensional state space and a 
discrete system with finite /^-dimensional Hilbert space. This 
situation is particularly relevant in quantum communication, 
where many advantages come from the use of interfaces be- 
tween light fields and matter- like systems [1], which are the 
basis of recent important experimental demonstrations UHl]. 
Ground-breaking architectures for communication, such as 
quantum repeaters [2], rely on interfaces and their efficient 
implementation is based on the availability of such techno- 
logical primitives. A "desideratum" of any reliable interface 
would be its ability to connect systems of different dimen- 
sionality so as to transfer important physical features (such as 
entanglement) across the interfaced systems while faithfully 
respecting their inherent structure. Among the proposals for 
quantum interfaces put forward so far US], those aiming at 
transferring entanglement from a light field encoding a CV 
state to a qubit system [4, 5] are appealing for implementa- 
tions in a number of physical settings, ranging from cavity- 
and circuit-quantum electrodynamics (QED) @, @] to polar 
molecules close to superconducting resonators |7|] or quantum 
dots and color-centers in diamonds in defect-microcavities 
and photonic crystals. Yet, a problem in dealing with mis- 
matched system dimensionality is the loss of any clear relation 
between the purity and entanglement properties of the CV re- 
source and those of an addressed two-qubit system. Control- 
ling that the state hierarchy properties of the fields are inher- 
ited by the qubits is thus difficult. 

In this paper we discuss a simple and flexible CV-to-qubit 
map, based on the use of handy Gaussian-light resources, 
which works as a powerful tool for quantum state-engineering 
of the steady- state of two remote qubits. The rich variety of 



mapped two-qubit states faithfully respects the state hierar- 
chy induced by the degrees of local/global purities and en- 
tanglement, whose structure is inherited from the Gaussian 
resources. Extremal bipartite Gaussian states (in terms of en- 
tanglement and mixedness) give rise to equally extremal two- 
qubit states. Our results thus shed new light on the mecha- 
nisms at work at the interface between discrete and CV sys- 
tems. Pragmatically, we show that our map represents an im- 
portant step in the long- sought task of distributing channels 
for quantum communication [1] and is implementable with 
resources available in many optical labs and relying on exist- 
ing mature technology in cavity- and circuit-QED laSl- 

The map - To keep the discussion as concrete as possi- 
ble, in the following we will adopt the language of cav- 
ity QED, although our scheme is completely general and 
independent of any specific physical setting. Two remote 
single-mode cavities contain a qubit each and are coupled 
to a broad-band two-mode driving field prepared in an en- 
tangled state £12 Qubit A is coupled to cavity 
mode 1 via the well-known Jaynes-Cummings Hamiltonian 
Ha \ -co{q\&\ + cr 2 A p\), where {cjj,Pj} are the quadratures of 
mode j = 1,2 and cr[ (2) is the x (y) Pauli matrix of qubit 
k=A, B. A similar model holds for qubit B and mode 2. Here, 
{|0), \ l)}k are the logic states of qubit k and oj is the Rabi fre- 
quency. The driven cavities leak photons at rate k and em- 
body local environments for the qubits. We are interested in 
the control and manipulation of the entanglement and purity 
properties of the two-qubit system by means of its interac- 
tion with the bosonic field. States of CV systems can be 
described in terms of a (generally infinite) hierarchy of mo- 
ments of their quadratures. Without loss of generality, the 
first moments are hereafter set to zero as they play no role 
in characterizing entanglement and mixedness. The second 
moments are used in the covariance matrix (CM) V12 of el- 
ements Tr[£i 2 {x a , xp\\ (a,j3 = 1,..,4) with x = {qi,pi, q^pi]. 
Gaussian states are completely specified by the knowledge 
of the first and second moments [9]. The CM of any two- 
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mode state can be brought in the form V12 = I qt y I w * m 

Vi = diag[#, a] (V2 = diag[Z?, b]) and C12 = diag[c+, c_]. In or- 
der for £12 to be a physical state, we must satisfy both V12 > 
and the uncertainty inequality Vi2+/H>0 with O = Q) k=1 2 ^ 
the symplectic matrix. In the bad-cavity limit where modes 1 
and 2 reach their stationary state sooner than any changes in 
the qubit- field subsystems and within the first Born-Markov 
approximation, the state of qubits A and B evolves according 
to the master equation (ME) [5] 



o t Qab = _T dj k (djQ AB d k -\d k 6j, q ab }/2), 
j,k=i 



where Oj=& 1 A 
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for 
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j=3,4. The Kossakowski matrix D of elements reads 
D=y(Vi2+/n), where y=2oj 2 /K is the effective qubit-field 
coupling strength, k is the cavity decay rate, and r-yt is the 
dimensionless time 

Eh. 

The map in Eq. (Q]) is completely 
positive iff D>0 01, which is equivalent to the uncertainty 
principle for the field CM. Thus the mapping holds for any 
(Gaussian or non-Gaussian) two-mode state with bona fide 
CM. Here, we consider only Gaussian states (generally mixed 
and asymmetric), bearing in mind that our results hold also for 
maps driven by non-Gaussian states with the same CM 11011 - 

In Ref. BI2I1 we study the dynamics of the two-qubit 
system. Here we deal with the features of £as(°°X 
which is found by setting d T g A B(r) = in the ME, 
in relation to q u . By calling Qij,ki=(ij\QAB(°°)\kl) 
(with \ij) and \kl) states of the two-qubit basis and 
i,j,k,l = 0,l), we get £00,00 = [(ab-a-b)z+(a + b) 2 ]/S, 
£01,01 = £00,00 + 2[<2z - (a + b) 2 ]/6, £10,10 
£00,00 + 2[bz - (a + b) 2 ]/S, £00,11 = 011,00 = 2(0 + b)(c--c+)/6, 
£01,10 =£10,01 =2(a + b)(c-+c+)/6 with d = 4abz, 

£11,11 = 1 -£00,00 -£01,01 -£10,10, andz = (a + Z?) 2 -2(c 2 +c 2 ). 
Degrees of entropy and entanglement - The mixedness (or 
lack of purity) of the state £ of a /^-dimensional system can be 
quantified by the linear entropy S = [D/(D- l)](l-Tr£ 2 ), 
ranging from (pure states) to 1 (totally mixed states). 
For a CV Gaussian state £12 with CM V12 one has 
Sl(£i2)=1-1/ VdetVi2. Both for two-qubit and two- 
mode Gaussian states, separability is equivalent to positivity 
of the partially transposed density matrix fl3ll . The degree 
of violation of such a criterion provides an entanglement 
monotone, the negativity N(qab ) = max{0, | \g T ^ B 1 1 1 - 1 } Il4ll , 
where T A stands for partial transposition with respect to 
qubit A and || • ||_ is the trace-norm. For a two-mode Gaus- 
sian state £12 one has 7V(£i2)=max{0, (l-v~ 2 )/v~ 2 }, where 
v7 2 is the smallest symplectic eigenvalue of the partially 
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(l/V2)[A 12 -(A 2 2 -4detVi 2 )^ 



transposed CM, 

with A 12 = detVi+detV 2 -2detCi2. 
sian state is entangled if and only if v~ 2 < 1. 



A two-mode Gaus- 
We 

s + d, b = s- d and 



parameterize V12 by setting a 

V(/rf -h d ) 2 -4g 2 ± VCA -h d f -4g 2 2 A 

c± = , 1 \ 2 with h d = (2d A +g)(A+\) and 

f x = 4x 2 + (g 2 + \)(A-\)l2 (x = d,s) UM- The Gaussian state 
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FIG. 1: (Color online), (a) S l (qab) versus S L (gu) for 20000 random 
states. Darker (lighter) dots denote entangled (separable) qubit states. 
Points below the dashed line show states purified by our map. (b) 
N{qab) versus normalized field negativity for 15000 random states. 



£12 (with purity Sl(£i2) = 1 ~ 8 l ) is physical and entangled 
for s>l,\d\<s-l,2\d\ + l<g<2s-l and 
where s and d determine the properties of the reduced states 
£i ( 2) having CM Vi (2 ) according to S L (g k )=l-[s+(-l) k ~ l d]~ l 
with (k - 1,2). At set values of g, d and s, N(gn) grows 
monotonically with X. 

What properties of the two-mode state are transferred to 
the two-qubit system? Marginal properties are faithfully re- 
produced as S L (g A (B)) = S L (£i(2))[2 - S L (gi i2) )], where g A{B ) 
is the reduced state of qubit A (B). The proportional- 
ity between Sl(qa(b)) and Sl(£i(2)) entails that the state- 
symmetry is preserved by the map. However, one finds that 
S l (Qab) = 1 - l/(3tf 2 ) - l/(3b 2 ) - I6(a + b) 2 f/(3S 2 ) (with 
^ - [( a + b) 2 +4(c 2 + + c 2 )]), which shows that the two-qubit 
mixedness is not a simple function of g alone, but depends 
nontrivially also on s, d, and X. In particular, in the al- 
lowed range of parameters, S l(Qab) increases with the global 
field mixedness (i.e. with g) and its mean energy (param- 
eterized by s), but decreases with X: larger input entangle- 
ment at given entropy results in qubit states of higher purity. 
Notwithstanding the effective non-unitary dynamics, when 
£12 is pure (which occurs when a-b and c+= - c-= Va 2 - 1) 
Qab(°°) is also pure. In general, for a given field mixed- 
ness, Sl(@ab) cannot vary unconstrained. Analytically, we 
find that S l(Qab) admits tight upper and lower bounds that 
depend on the field mixedness Sl(Qu)- This is further con- 
firmed by random numerical sampling [cfr. Fig. \Ha)]. The 
maximum of S l(Qab) at a given g is found by optimiz- 
ing over s, d and X. This implies taking X- -1 and 
d=0. The corresponding two-qubit states tend asymptotically 
to the Werner state g^ B = p\0~) AB (0~\ + (1 - p)t/4, where 
\0-) AB = (100) - |11>W V2 and p=2/(l + g 2 ), for which 
S^ x (qab)= 1-4[1 + (S l (0i 2 ) - l)" 2 ]" 2 . On the other hand, 
S™ m (g A B) at a given g is obtained for s = (g + l)/2, X=l 
and d- (g - I) 12. Such S™ n (QAB) is achieved by the product 
states [(g-l)|01X01|+(g+l)|HXH|]/(2g), for which we have 
S™ n (g A B) = (2 / 3) S l(qu)[2 - S l(q 12)] • The protocol can also 
be used to realize state purification: one finds many mapped 
Qab 's whose mixedness is smaller than the input S l(Q\i), even 
for totally mixed fields [see Fig.QJa)]. Noticeably, two-qubit 
states obtained from highly mixed £12 's are separable while 
entanglement arises in the region of moderate mixednesses. 
By introducing y/= 2(fl+ ^^- ) and p = z(a - b)(\ + t_ 2 ) 1/2 , the 
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two-qubit negativity at the steady- state reads 

N(g A B)=m<ix{0,(2/d)[(a + b) 2 -d/4+p]}. (2) 

In the physically allowed range of parameters, we have 
d g N(QAB)^Q, d s,aN(Qab) ^ while the dependence on d is 
non-monotonic. This implies that the two-qubit entangle- 
ment increases with the marginal entropies, decreases with the 
global mixedness and, at fixed global and marginal entropies, 
increases with A. These are the very same patterns followed 
by N(qu) and therefore the map in Eq. ([7]) fully preserves 
the qualitative structure of bipartite entanglement. Quantita- 
tively, however, N(qu) does not determine directly the two- 
qubit negativity. For a given value of N(qu), the correspond- 
ing Qab range from separable to highly entangled. The behav- 
ior of N(qab) versus the negativity of randomly sampled CMs 
V12 is reported in Fig. [Tib). While N(qab) can vanish for ar- 
bitrarily entangled Q\2 's, we find a maximum of the two-qubit 
entanglement N™ x (q ab ) = 1 -2[(1 + N(g u )) 2 + l]" 1 that is 
achieved by pure states (both for fields and for qubits). 
Entanglement versus global mixedness - Maximally entan- 
gled two-qubit mixed states (MEMS) are defined as those 
maximizing a given entanglement measure at any fixed value 
of the global mixedness [16]. In the {Sl(qab), N(qab)} space, 
MEMS include the family of Werner states. The correspond- 
ing minimally entangled mixed states are just separable states. 
The Gaussian counterparts to MEMS (GMEMS) are two- 
mode mixed states with infinite entanglement, such as the 
two-mode squeezed thermal states, with a-b- ygcosh(2r), 
c+ = -c- = y^sinh(2r) in the limit r —> oo. We remark that, 
given all possible CMs, our scheme does not generate ev- 
ery possible two-qubit state. However, the set of Qab's that 
can be engineered by our process does fill the entire region 
of the {Sl(@ab), N(qab)} diagram physically allowed to two- 
qubit states. This result is illustrated in Fig. [2a) where we 
report the diagnostics of Qab's obtained from random g u 's. 
The upper bound to the physically allowed region includes 
states gJ B for which Af max (£ Afi )=[-l+3 yj\ -S L (g AB )]/2 if 
S l(Qab)<%/9, and zero otherwise. A direct way to obtain 
such states is by maximizing A, minimizing d and setting 
g = 1 /p, which fixes the global purity of the two-mode re- 
source to be equal to the |®~> component of q^ b . Interest- 
ingly, the field state associated to such parameters is precisely 
a GMEMS with g=l/p. The mapped two-qubit state con- 
verges to the corresponding boundary state when the squeez- 
ing in Qu is large. We can thus engineer MEMS of tunable 
entanglement/purity by adjusting purity and squeezing in V12. 
In Ref. LI 2] we present an additional study of N(qab) against 
the marginal entropies of the resource. 

Entanglement versus global and marginal entropies.- A 
very refined characterization of entanglement is possible 
in the space of global and marginal entropies, where all 
the entangled two-mode Gaussian states lie in a narrow 
region bounded by Gaussian least-entangled and Gaussian 
maximally-entangled mixed states (GLEMS and GMEMS, re- 
spectively) [1]. These are achieved, at given s, d, and g (fix- 
ing the entropies) for A = -l (GLEMS) and A- 1 (GMEMS), 
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FIG. 2: (Color online), (a) N(qab) versus S L (g AB ) of two-qubit states 
obtained from 20000 random entangled two-mode Gaussian states 
Qi2 9 s. The upper boundary (MEMS) includes Werner states, (b) 
N(qab) versus S\ oc and S for states obtained using 10 4 random sym- 
metric entangled Q\ 2 's. (c) The same states as in (b) and the surfaces 
of maximum and minimum negativity at fixed entropies. 



respectively. Thus, by accessing only the restricted set of pa- 
rameters that determine the marginal and global purities of a 
two-mode Gaussian state, one can pin-down its entanglement 
and bound it with the corresponding GLEMS and GMEMS. 
Surprisingly, no exact two-qubit counterpart to this structure 
was known thus far. We will now show that, via Eq. (Q]), 
the two-qubit states qab inherit and enhance the properties 
of such CV states. As A e [-1, 1], the existence of GLEMS 
and GMEMS stems from observing that d^N{Q\2)\ s ,d, g ^ 0. 
A similar property holds for qab- For ease of notation, we 
write S=S l (qab) and N q =N{QAB)- From Eqs. © it fol- 
lows that N q is a function of s, d, A and g(s, d,A,S) such that 

dA^f q \s,d,s = d A N q \ s ,d, g + d g N q \ S 4,xdxg\s,d,s = d*Mq\s,d,g + 
dgN q \ S 4,A(d A S\ s ,d, g /d g S\ S 4 jA ) > 0. This shows the existence 
of two-qubit least-entangled and maximally-entangled mixed 
states (QLEMS and QMEMS, respectively), at fixed global 
and marginal entropies. These are obtained by mapping of 
GLEMS and GMEMS, respectively. The values of s and 
d are set by S^Qk) (k = A, B), while g is determined by the 
marginal and global entropies. For symmetric QMEMS with 
A=l, S l(qjc) = S\oc = 1-s~ 2 and global mixedness S , one has 
g = 3/[l-5ioc+ V4-95+5ioc(4+5ioc)]. The corresponding 
yymax ^ see ^ 2]) fixes the upper bound for all two-qubit 
states obtainable by our process and compatible with the given 
entropies. An analogous analysis holds for QLEMS. 

The behavior of the negativity versus global and marginal 
entropies shows that all the mapped two-qubit states lie in a 
quasi-bidimensional region [see Fig. [2b)]. Even for mixed 
states, Nab is almost perfectly a function of the global and 
marginal entropies alone. By superimposing the boundary 
curves corresponding to QLEMS and QMEMS to the numeri- 
cal analysis in Fig. Ob) we see that all the randomly generated 
states Qab's accumulate in the tight interval between maxi- 
mum and minimum negativity [cfr. Fig.Oc)]. Numerically, 
the negativities of QLEMS and QMEMS differ by less than 
0.04 e-bits close to the separability point, while in the region 
of larger entanglement they practically coincide. Therefore 
the field-qubit interface defines, within the set of entangled 
two-qubit states, the analogues to two-mode GMEMS and 
GLEMS. Moreover, the situation typical of the Gaussian case 
is enhanced in the case of qubits, since in the latter the gap be- 
tween maximal and minimal entanglement is even narrower. 
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The introduced mode-qubit dynamical interface preserves the 
hierarchy of entangled states in the entropic space: GMEMS 
(GLEMS) are mapped into QMEMS (QLEMS). Furthermore, 
qab weakly depends on X so that N(qab) can be accurately 
determined only controlling the engineered entropies of Qn. 
Practical considerations - Cavity QED is a natural setting 
for the implementation of the proposed scheme, due to the 
availability of a variety of CV resources LI 8]. The feeding 
of an optical cavity containing a trapped atom with squeezed 
light (bandwith of about 12 MHz) and the corresponding 
controlled light-atom interaction have been experimentally 
demonstrated [8] for (K,a))/2n - (70,20) MHz. This gives 
an effective coupling rate y ^ 10 MHz. The engineering 
of MEMS from approximated GMEMS produced at high 
squeezing is thus feasible. The swift progress in circuit-QED 
makes also such setting appealing [6] . A qubit is embodied by 
a superconducting quantum interference device (SQUID) at 
the charge-degeneracy point. The qubit transition-energy can 
be set by an in situ magnetic flux that modulates the Joseph- 
son energy so as to adjust the qubit- to-light coupling. Each 
qubit is integrated in a full/half- wave waveguide split by in- 
put/output capacitances: we can thus consider two discon- 
nected regions of a coplanar waveguide, joined via indepen- 
dent input/output capacitive lines for the injection/leakage of 
the field resource or via a large Josephson junction. The res- 
onator quality-factor is typically well within a range appro- 
priate to our scheme (10 2 to 10 6 with qj~0A GHz). For the 
frequencies involved in this setup, this would allow the real- 
ization of the regime studied here. Very recently, the ability 
to perform complete state tomography of two such qubits has 
been experimentally demonstrated 11191] . opening up the pos- 
sibility for preparation and characterization of the qubit state. 
Although the ability to experimentally engineer microwave- 
field states are inferior to their optical counterparts, squeezed- 
vacuum/thermal states have been produced by Josephson pa- 
rameteric amplifiers embodied by large junctions H20I1 . which 
is a very promising step towards the preparation of microwave 
states belonging to the classes studied here. Our scheme may 
be used to entangle collections of remote matter-like qubits, 
so as to achieve key resources in quantum technology. 
Conclusions- We have discussed a CV-to-qubit map that en- 
gineers two-qubit states spanning the region of physically al- 
lowed quantum correlations at fixed values of local and global 
entropies. It incorporates the most relevant sources of noise 
affecting a CV-to-qubit interface. Our results assure the re- 
alistic possibility for non-demanding production and control 
of qubit states using off-line preparation of CV entangled re- 
sources and linear local interactions. The relations found be- 
tween two-mode fields and two-qubit states make our scheme 
a basic predictive tool for light-matter entanglement transfer 
and related implementations in quantum technology. 
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Supplementary Material 



QUANTITATIVE ANALYSIS OF THE DYNAMICAL 
CV-TO-QUBIT MAP 

In this Section we provide further details on the map em- 
bodied by the master equation (ME) in Eq. (2) of the main 
paper. Rather than discussing the steady state properties of 
the mapped two-qubit state, here we address the full dynami- 
cal evolution of an initial preparation of qubits A and B. 

When written in the computational basis 
{|00>,|01>,|10>,|11>}ab, Qab can be partitioned as 
Qab(t) = Q x AB (j) + q° b (t), where 













Q X ab(t) = 


* * 

* * 


> Qab^T) = 


* * 

* * 


• (S-l) 













Here, the symbol * is used to denote any potentially non- 
zero matrix entry. Obviously, Q° AB (j) is not a density ma- 
trix. It is straightforward to check that the Liouvillian £ in 
Eq. (2) of the main paper keeps q\ b (t) disjoint from Q° b (t), 
so that starting from a state having elements only in £*#(t), 
we are sure that Q° b (t) = Vr. This is clearly the case for 
Qab(0) = |00)Afi(00|, and this observation yields a great simpli- 
fication enabling us to swiftly determine the dynamics of the 
entanglement transferred to the qubits. Moreover, such an ini- 
tial condition is the most favorable one to achieve the highest 
possible qubit entanglement in the model under investigation. 

Starting from the ME in Eq. (2) of the manuscript, it 
is straightforward to write down the B loch-like equations 
for the evolution of the two-qubit density matrix elements 
and solve them against the dimensionless interaction time 
r for any specific assignment of the resource's covariance 
matrix (CM). Fig. Oa) shows the typical r-dependent be- 
havior of such elements, where we have used the notation 
QijM = {ij\QAB(j)\kl) (with \ij) and \kl) states of the two-qubit 
basis and i, j, k, I- 0, 1), similarly to what has been done in the 
main paper. Similar trends are found for any allowed choice 
of parameters in V12. Interesting information can be extracted 
from the dynamical behavior of Qab(j). First, by fixing s, d 
and g in V12, one can easily see that the dynamical evolution 
of the negativity of the two-qubit system keeps N(qab) within 
a very narrow range of values, as A is varied (we remind that 
A e [-1,1]). This is in line with the behavior highlighted in 
Fig. 2(b) and (c) of the main paper, although here we are not 
imposing any restriction to the local/global entropies of the 
Gaussian resource state. Fig. Ob) shows, in this sense, an un- 
usual scenario where the difference between the negativity at 
A = -1 and the corresponding value at A = 1 isabout 30%. 
Typically the top-most curve would be less than 5% far from 
the bottom one. 

The information gathered throughout the dynamical evo- 
lution of the qubit density matrix allows one to infer how 



the two-qubit system approaches its steady state, for a given 
choice of the input CM. By using r as a sort of curvilinear 
abscissa of the negativity- versus-linear entropy functions, we 
can determine the trajectories of the two-qubit density ma- 
trices up to the boundary curve accommodating maximally 
entangled mixed states and Werner states [see Fig.lU. A sim- 
ilar study can be conducted in order to infer the dynamical 
mapped- state entanglement against the resource negativity (as 
in Fig. 0). We refer to the captions of such figures for further 
details on the simulations. 



ENTANGLEMENT VERSUS MARGINAL MIXEDNESSES. 

Here we classify the bipartite entanglement in the marginal- 
entropy space. Such a characterization is exact for pure states. 
In general, bounds on mixed-state entanglement can be de- 
rived for all mixed states compatible with given marginals. 
For two qubits, the upper bound to the physical set of 
states in this diagram defines maximally entangled states at 
fixed marginal mixednesses (MEMMS) fcj], while separable 
(product) states exist for any pair of marginals [3]. For 
equal marginal entropies, MEMMS are pure states. Their 
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FIG. 3: (Color online), (a) Behavior of the density matrix elements 
against the dimensionless interaction time r for (randomly taken) s = 
1.774, d = 0.07, g = 1.448,7 = 0.1 and A = 1. (b) Dynamical 
behavior of the negativity N(qab) against r for the same values as in 
panel (a) but A, which is taken to grow from -1 to 1 in steps of 0.4. 
The lower-most (top-most) curve is for A = -1 (A = 1). The dashed 
horizontal lines show the steady- state values. 
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FIG. 4: (Color online). Typical trajectories of the mapped two-qubit 
states g AB in the {S {qab), N l {q A b)} plane. The dashed boundary curve 
embodies maximally entangled mixed states (MEMS) and Werner 
states q™ b . The dots show the evolution of negativity and global lin- 
ear entropy of the two-qubit state for five different choices of the 
Gaussian resource parameters s, g, d and A. The dimensionless in- 
teraction time r grows as indicated by the arrow. The final dot in 
each trajectory indicates the corresponding steady state. The trajec- 
tory superimposed to the MEMS curve is for a configuration of the 
parameters entering V12 which guarantees the asymptotic mapping 
to a g^ B state. 

Gaussian counterparts have been introduced in Ref. U] and 
dubbed Gaussian-MEMMS, or GMEMMS. Their CM is char- 
acterized by having g = 2\d\ + l. Correspondingly, any de- 
pendence of the CM on A disappears and one simply has 
c+ = ± V(l + max{<2, b})(- 1 + minja, b}) QD. 

A random- state investigation of the performances of our 
map in such local-entropy space shows that the mapped two- 
qubit steady states do not fill the whole region allowed to any 
two-qubit state with given marginals. In particular, any at- 
tempt to reproduce MEMMS results in unphysical parame- 




Af(f?i 2 )/(l+A/"(f?i2)) 

FIG. 5: (Color online). Typical trajectories of the mapped two-qubit 
states q ab in the plane spanned by the (normalized) negativity of 
the Gaussian resource N(qu) and the target qubit- state negativity 
N{qab)- The boundary curve follows the functional form found in 
the main paper. The dots along each vertical line show the evolution 
of the two-qubit negativity for random choices of the Gaussian re- 
source parameters s, g, d and A. The dimensionless interaction time 
r grows as indicated by the arrow (at r = 0, obviously, each two- 
qubit state is separable). The final dot in each trajectory indicates the 
corresponding steady state. 
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FIG. 6: (Color online). Two-qubit negativity against local marginal 
entropies. Physical two-qubit states lie below the top tent-like bound- 
ary (MEMMS). Such states are not reproducible using Gaussian 
states as these produce Qab's lying in the inner region, whose upper 
boundary contains states engineered via GMEMMS. 



ters for the driving field. In Fig. [6] we show such a numeri- 
cal exploration of the negativity N(qab) against the marginal 
entropies S l(Qa) and Sl(@b)- The attainable qubit state fill 
a restricted, pyramid-like region: their entanglement is never 
maximal at given marginals, except in the pure- state case of 
Sl(Qa) = Sl(qb). We now aim at characterizing analytically 
the upper boundary of this set. Recalling that the negativity is 
an increasing function of A and a decreasing function of g (at 
fixed a and b), we can conclude that the boundary is obtained 
by minimizing g, which gives g = 2\d\ + 1. Such values of the 
parameters are exactly those characterizing a GMEMMS. Al- 
though the full range of entanglement is not achievable in the 
space of the marginal mixedness, our map is such that the two- 
qubit states endowed with the maximum achievable entangle- 
ment at fixed marginals are those obtained by using GMEMS. 
Therefore, despite being unable to reproduce MEMMS, the 
maximum negativity achieved by our scheme corresponds to 
the images of GMEMMS. Fig. [6] provides further details on 
this point. We believe this is yet another clear indication of 
the powerful and faithful mapping embodied by by our sim- 
ple bilocal linear interaction model. 



ANALYTIC EXPRESSION FOR THE NEGATIVITY OF 
QMEMS 

We consider the space of negativity against global and 
marginal mixedness. The aim of this Section is to give an 
explicit form to the negativity of two-qubit most-entangled 
mixed states (QMEMS). By taking the case of a resource 
embodied by a Gaussian most entangled mixed states (or 
GMEMS) with A = 1 and imposing the symmetry conditions 
given by S L (g A ) = S L (g B ) = 1 -1/s (we call S L (g AB )=S the 
global linear entropy), it is straightforward to check that the 
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negativity of the extremal two-qubit states in this space reads 



jV max = I 
q 6 



-(2 + S loc ) + J(2 + S loc Y-9S 



+2 [l + 8S loc - Sl c -9S + (Sioc-1) ^(2 + S l0C f-9S 

(S-2) 

if 9S + 4(-2 + S loc )S loc < 4Sioc ^/l-^ioc + ^L and = 
otherwise. More cumbersome expressions are found for g 
and the corresponding minimal value N™ m valid for two-qubit 
least entangled mixed states (QLEMS). 
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